Quantum Hall plateau transition in graphene 
with spatially correlated random hopping 
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We investigate how the criticality of the quantum Hall plateau transition in disordered graphene 
differs from those in the ordinary quantum Hall systems, based on the honeycomb lattice with 
ripples modeled as random hoppings. The criticality of the graphene-specific n = Q Landau level 
is found to change dramatically to an anomalous, almost exact fixed point as soon as we make the 
random hopping spatially correlated over a few bond lengths. We attribute this to the preserved 
chiral symmetry and suppressed scattering between K and K' points in the Brillouin zone. The 
results suggest that a fixed point for random Dirac fermions with chiral symmetry can be realized 
in free-standing, clean graphene with ripples. 
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After the seminal observation of the anomalous quan- 
tum Hall effect (QHE) in graphene, [1 13, Q fascination 
expands with the graphene QHE. One crucial question 
that is not fully explored is: what exactly is the role of 
the chiral symmetry in the problem? This has to do with 
a most significant feature of double Dirac cones (at K and 
K' in the Brillouin zone) in graphene. Although a single 
Dirac cone would already imply a characteristic Landau 
level structure with the zero-energy level, if we really 
want to look at the effect of disorder on the graphene 
Landau levels, we have to go back to the honeycomb lat- 
tice for which the chiral (A-B sub-lattice) symmetry^, Q 
and the associated valley (K and K') degrees of freedom 
enter as an essential ingredient. The effect of disorder 
should then be sensitive to the nature of disorder, i.e., 
bond disorder or potential disorder, which determines 
the presence or otherwise of the chiral symmetryj^, 0,11], 
and whether the disorder is short-ranged or long-ranged, 
which controls the scattering between K and K' points. 

For Dirac fermions, effects of random gauge fields in- 
duced by ripples in the two-dimensional plane have been 
discussed[^ [13, El, and the stability of zero modes has 
been argued in terms of the index theorem and the chiral 
symmetry [l|, 0, [E Ell- More recently, the plateau-to- 
plateau transition for random Dirac fermions has been 
discussed, where the particle- hole symmetry is shown to 
make the zero-energy Landau level robust [13| . As for the 
criticality, however, the result shows nothing special 
about the n — Q Landau level, but this is obtained for 
a model of the Dirac fermions for which the randomness 
is introduced as a scalar random potential, so the chiral 
symmetry is degraded. 

On the other hand, the actual randomness in graphene, 
even when atomically clean, is known to have ripples, 
i.e, long-ranged corrugation of the graphene plane 9]. In 
fact, while a monolayer graphene naively contradicts with 
the well-known theorem that two-dimensional crystals 



should be thermodynamically unstable, one explanation 
attributes the stability to the ripples[14]. In this sense, 
we can take the disorder coming from ripples in graphene 
as an intrinsic disorder. Since the ripples consist of ran- 
dom bending of the honeycomb lattice, its main effect 
should be, in the tight-binding model, a modification of 
the hopping integral between neighboring sites Jii]. Thus 
the question amounts to: how does the QHE criticality 
behave for a model with random hopping on the hon- 
eycomb lattice. The random hopping is of fundamental 
theoretical interest as well, since a bond randomness pre- 
serves the chiral symmetry, so its effect, particularly on 
the criticality, is of crucial interest. The chiral symme- 
try indeed plays a fundamental role in graphene 0, [HI, 
which protects the gapless double Dirac cones as well as 
the existence of characteristic zero- modes with/without 
magnetic fields. A special importance of the chiral sym- 
metry in localization physics has also been discussed with 
a viewpoint of the universality [15i, ilQ, il7i] . 

Now, the length scale over which the lattice is warped 
should be reflected as the spatial correlation in the ran- 
dom hopping. In the case of the ordinary QHE systems, 
the importance of the spatial correlation of randomness 
has been discussed in various contexts, among which are 
the pioneering work by Ando and Uemura[18|, the levi- 
tation of the critical states in the lattice mo de l [Til . [20| , 
a multifractal analysis of critical wave functions [2 ll|. and 
plateau transitions in narrow wires (22]. In graphene, on 
the other hand, the range of disorder plays an unusually 
important role, since the range dominates the inter-valley 
(K-K') scattering. Hence we conceive here that it is im- 
perative to examine the honeycomb lattice (rather than 
an effective Dirac model) with bond randomness of var- 
ied correlation lengths in understanding the Hall plateau 
transition in graphene. 

This is exactly our motivation here to explore how the 
criticality in disordered graphene QHE transition, espe- 
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cially for n = Landau level, depends on (i) the symme- 
try and (ii) the range of randomness in the honeycomb 
lattice. As for the Hall conductivity which has a topolog- 
ical ori gin and mathematically a Chern number in units 
of e^/ft. 24, [2^ [2^, an unusually accurate and efhcient 
method is required for examining QHE around the Dirac 
point (band center) for random systems. Here we have 
adopted a non-Abelian extension of the Chern-number 
formalism as combined with a lattice-gauge technique. 
We shall show that, while the plateau transition for the 
n = Landau level has an ordinary critical behavior for 
the Mworrelated random bonds, the criticality changes 
dramatically to an anomalous, almost exact fixed point 
with a step-function-like plateau transition and a con- 
comitant delta-function-like Landau level, as soon as the 
spatial correlation in the random bonds exceeds only a 
few bond lengths. This can indeed be attributed to the 
preserved chiral symmetry, which is confirmed by adding 
site randomness to modify the symmetry. 

The tight-binding Hamiltonian for the honeycomb lat- 
tice is H = J tijC^^'-' c\cj, in standard notations, where 
the Peierls phase } is determined such that the sum 
of the phases around a hexagon is equal to the mag- 
netic flux —2'K(j} piercing the hexagon in units of the flux 
quantum = h/e. The spin degrees of freedom are 
neglected for simplicity. We introduce randomness in 
the nearest-neighbor transfer energy as tij — t -\- 6tij, 
where the disordered component dtij is assumed to be 
Gaussian distributed, P{St) = e'^^^ 1'^"''' j ^phw^ , with 
a variance o. Next we specify the spatial correlation 
r\ in the random components by requiring (StijStki) = 
((5t^)e~l'"'^^'''''' , where denotes the position of 
the bond tij, and () the ensemble average[22|. We take 
the X- and j/-axes as shown in FigH] for x Ly rectan- 
gular systems, where a typical spatial landscape of the 
random hopping is displayed. All lengths are measured 
hereafter in units of the bond length a in the honeycomb 
lattice. 

Even with such random transfers, the Hamiltonian re- 
spects the chiral symmetry, that is, there exists a local 
unitary operator 7 (with 7^ = 1), which anti-commutes 
with the Hamiltonian, {_ff, 7} = 0. In real space we can 
decompose the honeycomb lattice into two sub-lattices A 
and B, for which the fermion operators are transformed 
as 7Ci7~^ = sCi with s = -|-1(— 1) for i G A{B). Obvi- 
ously, this symmetry is destroyed by a potential disorder, 
while the random hopping preserves it, even in magnetic 
fields. Since the eigenstates appear in chiral pairs {ip, 7V' 
with eigenenergies ±iJ), it is clear that the zero-energy 
states are special. If a zero-energy state ip is not an eigen- 
state of the chiral operator, we can use the zero-energy 
chiral pairs, tp and jip, to make them eigenstates of the 
chiral operator, "fip± = i4'± with 'ip± — ipi^tp. Then all 
the zero modes are eigenstates of 7 with amplitudes re- 
siding only on one of the A and B sub-lattices. Hence the 




FIG. 1: (Color online) An example of the spatial landscape 
of the random components, St/a, in the hopping for a spatial 
correlation length rj/a = 5. 



topologically protected zero-energy Landau levels, partic- 
ularly their criticality, can be very sensitive to whether 
the disorder respects the chiral symmetry or not. 

Let us first look at the result for the density of states 
around E — 0. A key interest is how the n = 
Landau level is broadened by randomness as compared 
with n ^ levels. In the calculation of the density of 
states we adopt the Landau gauge for the correspond- 
ing bricklayer lattice [23] with periodic boundaries in y- 
dircction and armchair edges in x to remove the contri- 
bution from zigzag edges. The density of states (pi) = 
— '^■lmGii{E -\- vy)/NTT, is obtained in terms of the 
Green's function[27|, Gii{E + ie) = {i\{E - H + ie)-'^\i), 
where N is the total number of sites, and e a small 
imaginary part in energy to evaluate the Green func- 
tion numerically. We have performed the calculation for 
6.3 X 10"'* < e/t < 1.0 X 10-2 and confirmed that the 
anomaly at _E = described below is not affected by the 
value of e. 

The result for the density of states with the disorder 
strength a/t = 0.12 and a magnetic field 0/(/>o = 1/50 
for various values of the correlation length rj is shown in 
Fig. [21 It has been shown [J] that the n = and several 
adjacent Landau levels characteristic to the relativistic 
electrons are captured even with this value of (j), which, 
when directly translated, corresponds to a large magnetic 
field, so the model should be adequate for the analysis of 
the criticality at the n = Landau level. We can imme- 
diately see that the n = Landau level is anomalously 
sharp, but that the sharpness depends sensitively on the 
correlation length rj of the random hopping. More pre- 
cisely, as soon as we have ry/a > 1, the n = Landau level 
becomes remarkably sharp, while this does not occur for 
n 7^ Landau levels. Indeed, the shape of the n = 
Landau level for 77/a > 3 is delta-function-like within the 
numerical accuracy in that its shape coincides almost ex- 
actly with the Lorentzian density of states in the clean 
limit, p{E) = ^ E'-'+e'^ • -'^^ ^^^^ sense the n — Landau 
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FIG. 2: (Color online) The density of states for the spatially 
correlated random bonds for various values of the correlation 
length < 77/a < 2.0 with the strength of disorder fixed at 
(j/f = 0.12 in a magnetic field tp/(j>o = 1/50. The system size is 
L^/{VSa/2) = 5000, Ly/(Sa/2) = 100, and e/t = 6.25 x 10"''. 
Inset: The density of states for various values of the disorder 
strength a with a fixed ri/a — 3. 



level in the presence of the correlated bond randomness is 
delta-function- like for rj/a> 3. If we examine the depen- 
dence of the density of states on the disorder strength, 
with a fixed 77/a = 3 (Fig. 2, inset), we can confirm that 
this anomaly aX n = Q[E — Q) remains insensitive to 
the disorder strength as far as 77/a > 1, whereas other 
Landau levels are broadened by disorder. 

We now turn to the quantized Hall plateau transition. 
The Hall conductivity Uxy is related to the Chern num- 
ber nc as Gxy = nc(e^//i) provided that an energy gap 



exists above the Fermi energy [2J, |25|, |26[. In random 
systems the Chern number differs from sample to sam- 
ple, so we should look at the ensemble-averaged quantity 
for each energy bin, which gives the Hall conductance as 
a function of £; [ii,!!^. 

Since the Hall current is dissipationless, not only the 
state near the Fermi energy but all the filled states con- 
tribute. A speciality of the graphene QHE is that the 
region of interest is around E = 0, which implies that we 
have to question many Landau levels below the Fermi en- 
ergy whose contributions almost cancel with each other 
to a value of order unity. So we obviously confront a nu- 
merically difficult situation, especially if we want to look 
at a criticality around E = 0. We have previously shown 
that such a situation can be treated with a non-Abelian 
formulation of the Hall conductivity as a Chern number 
for multi-dimensional multiplets of fermions0, H^. The 
Berry connection is then defined as a matrix, which is 
spanned by the Landau sub-bands, split by the random- 
ness. We can then adopt an extended unit cell for each 
realization of the randomness to apply the formula 25 1. 
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FIG. 3: (Color online) The Hall conductivity (Chern number) 
against the Fermi energy E/t for spatially correlated random 
bonds for the correlation length 77/0 = 1.5(a) and 77/a = 0(b). 
We have a disorder strength a jt = 0.12, a magnetic field 
(?!)/(/>o ~ 1/50, a system size Lj:/{\/ia/2) = Ly/{'ia/2) = 20, 
= 10 and an average over 300 samples. 



The required stable energy gap at the Fermi energy is 
mostly guaranteed by the level repulsion in finite, ran- 
dom systems. Level crossings below the Fermi energy 
do not cause any problem either in this formulation. 
In the numerical evaluation of the topological numbers, 
a technique developed in the lattice gauge theory has 
turned out to be useful, which is a two-dimensional 
generalization of the King-Smith- Vanderbilt formula for 
polarization 3^, 3l| . For this we employ twisted bound- 
ary conditions, il){x + Lx,y) = e^'^'"ip{x,y),tfj{x,y + 
Ly) = e"t'y^p{x,y), where <j>x{y) = 27rnx(y)/N with 
"■a:(y) — 0, 1, 2, . . . , — 1 being discretized phases. 
Here the string gauge is used to treat weak niag- 
netic fields with the twisted boundary condition 1 1911 . 
The Chern number is then evaluated as [3, 0, l3Cll | 
«c - 2^ E<^ arg(det UiU^+^^- [uf+'^^- Uf]*), where 
^x{v) ~ [*(^ + '^4>x{y))]^'^{4>) with a set of eigenstates 
*(</>) — {\ipi{cf))), • • • , \ipMi4>))) below the Fermi energy 
and A</)^ = (27r/A^,0), A(f>y ^ {0,2t:/N). 

The Chern number averaged over 300 realizations of 
randomness for a bond disorder strength a/t = 0.12 is 
shown in Figs. [3] and [?] plotted within the energy region 
{\E/t\ < 1), where the plateaus have the Dirac behavior, 
{2m + l)(e^//i) with m an integer^. K we first look at 
Fig. [3] for the correlated disorder with ?7/a — 1.5, we 
immediately notice that the plateau transition between 
nc = —1 and 1 around E = is anomalously abrupt, i.e., 
nc behaves like a step function, in sharp contrast with 
other transitions for n ^ 0. For the uncorrelated bond 
randomness 77/a = 0, on the other hand, the transition 
for 71 = is as smeared as those for 7i 7^ 0. The anoma- 
lously sharp step for ti = agrees with the anomalously 
sharp 71 = Landau level seen in the density of states. 
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FIG. 4: (Color online) The Hall conductivity (Chern num- 
ber) as a function of the Fermi energy E/t for spatially cor- 
related random bonds with r^/a = 1.5 for two system sizes, 
L^/{VSa/2) = Ly/{3a/2) = 20(a) and = 10(b). We have a 
disorder strength a/t = 0.12, a magnetic field (j}/(j}o ~ 1/50 
with an average over 300 samples. Insets: Plateau tran- 
sitions for n = (c) and n — —1 (d) when we add a 
potential disorder introduced as random site-energies uni- 
formly distributed over [—w/2,w/2] and spatially uncorre- 
lated, and the curves are for w/t — 0,0.05,0.1,0.2,0.3,0.4 
with L^/{V3a/2) = Ly/{3a/2) = 10. 



We can confirm the behavior for n = is indeed un- 
usual by looking at a system-size dependence for a cor- 
related disorder in Fig. [U in which the results for two 
sizes coincide with each other within the numerical accu- 
racy for the transition at E — 0, while other transitions 
exhibit the usual behavior of narrower transition widths 
for larger systems. We can further test our picture that 
the anomalous behavior is connected to the preserved 
chiral symmetry. For this purpose, we have added a po- 
tential disorder to the bond disorder (Fig. 4, insets). The 
result clearly shows that the addition of a potential disor- 
der that destroys the chiral symmetry does wash out the 
anomalous (step-function-like) transition at E' = into 
a normal behavior (Fig. 4(c)), while other transitions re- 
main essentially the same (Fig. 4(d)). 

In summary we have revealed that the quantum Hall 
transition at i? = is anomalously sensitive to the spa- 
tial correlation of the random bonds, where concomi- 
tantly with the Landau level width, it becomes exact 
fixed-point-like as soon as the correlation length exceeds 
a few times the bond length. This sharply contrasts with 
the case of the generic random Dirac fermions, where the 
broadening of the Landau level occurs also for E = 0. 
The singular behavior may correspond to the fixed point 
for the random Dirac fermions with chiral symmetry dis- 
cussed by Ludwig et al[15|, where a generic instability of 
the fixed point is discussed. 

Experimentally, the length scale of ripples is estimated 
to be several nanometers 0, Since this is much 



greater than the correlation length adopted here, the 
bond disorder by such ripples should not broaden the 
n = graphene Landau level. Conversely, the broad- 
ening at the n = level as observed in experiments 
should be caused by other types of disorder, such as po- 
tential disorders by charged impurities [lo[ . The message 
here amounts to that the fixed-point behavior should be 
experimentally observed in free-standing clean graphene 
samples where the ripple is the only disorder. 
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